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SUMMARY 


A study of the motion of a personal-owner or liaison airplane 
through the stall and into the incipient spin has been made by cata- 
pulting a dynamic model into still air. The results showed that the 
model became unstalled in inverted flight after its initial rotation of 
substantially l80° in roll. Since the rates of rotation were low, it 
may be inferred that termination of the incipient spin would be most 
readily obtained by proper control movements at this point. 


INTRODUCTION 


Statistics indicate that a large percentage of all fatal private 
flying accidents in the past have occurred because of pilot error; in 
more than half of these accidents, the airplane stall has been involved 
(ref. l) . Oftentimes, accidents of this nature are referred to as 
"stall-spin" accidents and they occur primarily in the incipient phase 
of a spin, that is, in that portion of the motion immediately following 
the stall and before a spin can fully develop. In addition, spin demon- 
strations of some airplanes for the armed services require recoveries 
from only two -turn spins which also may not be fully developed motions . 
Data on the incipient spin are not available although much effort has 
been expended in the Langley 20-foot free- spinning tunnel toward an 
understanding of the developed spin. The problems associated with the 
incipient phase of the spin may be quite different from those associated 
with the fully developed spin and recovery therefrom. 

A technique is being developed with free dynamic models, similar 
to those used in spin testing, to study the incipient spin. The model 
is catapulted into still air in such a manner that stalled flight occurs. 
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and the ensuing motion at and beyond the stall is studied by use of 
motion-picture records from which the angles of attack and sideslip and 
the linear and angular velocities are determined. Scale effect could 
greatly influence the motion at and just beyond the stall although it 
has been shown to be of only secondary importance in developed spins. 

The technique described, therefore, may be used primarily to isolate the 
effects of various parameters rather than to indicate specific incipient- 
spin characteristics. The initial investigation was performed to deter- 
mine the motion in the incipient spin of a model of a low-wing, single- 
vertical-tail airplane typical of a present-day four-place personal-owner 
or liaison airplane and to determine the extent to which the motion may 
repeat itself. The results of this investigation are reported herein. 


SYMBOLS 


The analysis of the motion was made with respect to the body system 
of axes . A diagram of these axes showing various angles pertinent to 
the analysis is given in figure 1. 


X,Y,Z 


earth axes, a system of mutually perpendicular axes 
through airplane center of gravity always per- 
pendicular and parallel to earth's surface and of 
constant direction 


x B> y B> z B 

x r R > x rT ,X r C g 


yr R' yr T ,yr cg 


z r£> z r T > z r 


eg 


X R- ,y R ,Z R 

x 'juy , p> z p 


body axes 

distance along X earth axis from reference point to 

right wing tip, tail, and center of gravity, respectively 

distance along Y earth axis from reference point to 

right wing tip, tail, and center of gravity, respectively 

distance along Z earth axis from reference point to 

right wing tip, tail, and center of gravity, respectively 

distance of right wing tip from center of gravity along 
X, Y, and Z earth axes, respectively 

distance of tail from center of gravity along X, Y, 
and Z earth axes, respectively 


b 

In 


span of model, ft 

tail length (distance between center of gravity and inter- 
section of Xg with trailing edge of rudder) , ft 
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c 

mean aerodynamic chord, ft 

x/c 

ratio of distance of center of gravity rearward of leading 
edge of mean aerodynamic chord to mean aerodynamic chord 

z/c 

ratio of distance between center of gravity and fuselage 
reference line to mean aerodynamic chord (positive when 
center of gravity is below fuselage reference line) 

m 

mass of airplane, slugs 


moments of inertia about X B , Y B , and Zg axes, 
respectively, slug-feet 2 

\ - % 
mb 2 

inertia yawing-moment parameter 

' % 

mb 2 

inertia rolling-moment parameter 

Iz B " Ix B 
mb 2 

inertia pitching-moment parameter 

P 

air density, slugs/cu ft 

V 

resultant velocity of center of gravity, ft/sec 

u,v,v 

components of velocity V along X, Y, and Z earth 
axes, respectively, ft/sec 

U B' V B> W B 

components of velocity V along X B , Y B , and Zg axes, 
respectively, ft/ sec 

a 

w 

angle of attack, tan“^- — , deg 

U B 

P 

_ V-r, 

angle of sideslip, sin -1 deg 

h 

altitude loss, ft 

P 

rolling angular velocity about X fi , radians/sec 
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pitching angular velocity about Y B , radians/ sec 


r 

yawing angular velocity about Z B , radians/ sec 

OA and OC 

projections of X B and Y B , respectively, on 
earth X,Y plane (fig. l) 

OB 

position X B would take if airplane were rotated about 
Y B until X B lay in earth X,Y plane (fig. l) 

OD 

position Y b would take if airplane were rotated about 
X B until Y B lay in earth X,Y plane (fig. l) 

OE 

position Z B would take if airplane were rotated about 
Y b until X B lay in earth X,Y plane (fig. l) 

OF 

position Z B would take if airplane were rotated about 
X B until Y b lay in earth X,Y plane (fig. l) 

0 

angle that X B axis makes with earth X,Y plane, 
positive when nose is above plane (Euler's angle 0) 

0 

angle between Y B and line OD (fig. l) (Euler's 
angle 0) 

♦ 

angle between X earth axis and OA measured as shown 
in figure 1 from 0° to 36O 0 (Euler's angle \|r) 

♦ 

rate of change of \|r with time, radians/sec 

0 

rate of change of 0 with time, radians/sec 

0 

rate of change of 0 with time, radians/sec 


angle between Y earth axis and OC measured as shown 
in figure 1 from 0° to 36O 0 

0’ 

angle that Y B axis makes with earth X,Y plane posi- 
tive when right wing is below plane 

9 ' 

angle between X B and line OB (fig. l) 

♦ 

angle in X,Y plane between OA and OB or between 1 
and OD (fig. l) 
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APPARATUS AND METHODS 
Model 


The model used in this investigation was a configuration similar 
with respect to mass and dimensional characteristics to one investigated 
in reference 2 for fully developed spin and recovery characteristics. 

Based on the model size and the average dimensions obtained for a large 
number of personal-owner airplane designs, the model may be considered 
to be a 1/12.4 scale of a personal-owner or liaison airplane. The weight 
and moments of inertia of the model represent those of an airplane typical 
of this type. A three-view drawing of the model is shown in figure 2 and 
the mass and dimensional characteristics of a corresponding full-scale 
airplane are tabulated in table I. 

The model was constructed principally of balsa and reinforced with 
spruce and cedar. A photograph of the model is shown in figure 3 . In 
accordance with spin- tunnel practice, the propeller and landing gear were 
not simulated on the model. The model wing had an NACA 43012 airfoil 
section which at low Reynolds number has a sharp break in the lift curve 
at the stall. (See ref. 3*) 


Apparatus 


The tests were performed inside a building approximately 70 feet 
square and 60 feet high. The launching apparatus, shown in figure 4, 
was located 55 feet above the floor near one wall of the building and 
consisted primarily of an elastic chord which catapulted the model from 
its launching platform along a short track. The launching platform 
could be adjusted to set the model attitude and thus the launching angle 
of attack. An electronic timer was used to measure the model velocity 
just after it was launched. Fixed motion-picture cameras were located 
along the walls parallel and perpendicular to the initial launching 
direction to record the model motion. The cameras had overlapping fields 
so that the model was in view of at least two of the cameras, located at 
right angles to one another, throughout its flight. The motion-picture 
records were synchronized by use of a common timing light reflected into 
the camera fields. A large net for model retrieving purposes was hung 
from the wall opposite the launching apparatus. 


Testing Technique 

The model was catapulted from the launching apparatus at a speed 
slightly in excess of the stalling speed and at an angle of attack just 
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below the stall angle of attack. The elevator control of the model was 
set so as to pitch the model from the launching angle of attack up through 
the stall and the rudder was set to initiate a yawing motion and thus to 
precipitate a roll-off after stalling. 

As previously mentioned, photographic time histories of the model 
motion were made with a system of motion-picture cameras arranged so as 
to have the model in the field of two cameras with perpendicular focal 
axes at all times. When at an instant of time the model is in the field 
of two cameras located at right angles to one another, the space coordi- 
nates of the model can he determined. By this means, the model could he 
located at any instant from the time of catapult until the time the model 
struck the retrieving net. From a knowledge of the space coordinates 
and their instantaneous time values, the velocities of the model were 
obtained along the earth axes. Then, the attitude of the model at any 
instant being known, angle of attack, angle of sideslip, and linear and 
angular velocities were calculated. The derivation of equations for these 
calculations is presented in the appendix. 


Test Conditions and Accuracy 

For the tests, the model was ballasted with lead weights to repre- 
sent a corresponding airplane at an altitude of 5,000 feet (p = 0.0020L9 
slug/cu ft) . 

For this investigation, the controls of the model were set so as to 
cause the model to stall and enter a spin soon after launching. The 
settings of the controls were as follows: 

Elevator full up, deg 30 

Ailerons, deg 0 

Rudder fully deflected (in direction desired for roll-off), deg . • • 30 

The stalling angle and stalling speed of the model were determined 
from force tests of the model conducted on a strain-gage balance in the 
Langley 20-foot free-spinning tunnel. The Reynolds number of the present 
investigation was approximately 1 X 10 5 based on the mean aerodynamic 
chord of the model. 

Errors in the measured motion were magnified because graphical dif- 
ferentiation was required in reducing the data. For the current tests, 
the accuracy was estimated to be within the following limits: 
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a, deg *5 

0 , deg ±4 

u, v, w (full scale), ft/sec ± 1.3 

u-g, V B> sca l e )> ft/sec ±4.9 

p, q, r (full scale), radians/sec ±0.15 

Altitude loss (full scale), ft ±1 

Weight, percent _ ±1 

Center-of -gravity position, percent c ±1 

Moments of inertia, percent ±5 

Control settings, deg ±1 


RESULTS AND DISCUSSION 


The results of the tests are presented in figures 5 and 6 . The 
model test data have been converted to full-scale values for a corre- 
sponding airplane at a test altitude of 5,000 feet based on the premise 
that the model scale was l/l2.4. Film strips showing the motion of the 
model as it stalls and enters a spin are presented in figure 7 » 

The results as the model rolled to the right and entered an incipient 
spin are presented in figure 5 for an initial and a repeat test. Varia- 
tions with time of the angle of attack, angle of sideslip, altitude loss, 
and the rates of roll, pitch, and yaw are given in figure 5. Sketches of 
the approximate attitude of the model at various phases of the motion are 
also shown in figure 5. As was previously noted, the elevator was full 
up ( 3 O 0 ) and the rudder was full right ( 3 O 0 ) for these tests. Soon after 
the model stalled, it again became unstalled (time, approx. 1.75 sec), 
went inverted for a part of its motion, and assumed a very low angle of 
attack momentarily before again stalling and continuing toward a developed 
spin motion, such as reported in reference 2. After the initial stall and 
immediately after the model became unstalled, the rates of yaw and pitch 
were relatively small and the rate of roll began to decrease. There also 
was little loss of altitude up to this time and the results indicated that, 
even though the airplane may be inverted, a time soon after the roll-off 
has started appears to be a desirable time to attempt to terminate the 
motion. This conclusion is in agreement with pilots' experience that 
recoveries attempted soon after an airplane stalls and rolls off are 
more readily obtainable than those attempted from developed spins. 

The motion obtained may be qualitatively analyzed as follows: After 

being launched, the model immediately pitched up and yawed right because 
of the initial settings of the elevator and rudder, respectively. As the 
motion progressed, the sideslip angle (3 became negative because of the 
yawing motion; this sideslip angle tended to cause the model to roll off 
as a result of the dihedral effect. As was previously noted, the wing 
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had an NACA 4 3012 airfoil section for which the lift curve has a sharp 
break at the stall (ref. 3)} this effect undoubtedly augmented the model 
roll-off because of a loss in damping in roll beyond the stall. When 
the model again became unstalled, the damping in roll led to decreased 
rate of roll (p) until the stall reoccurred with the resulting increase 
in rate of roll. Changes in sideslip angle occurred as the model rolled 
and increments of positive and negative velocity components along the 
YB axis were obtained as the right or left wing, respectively, was 
inclined downward. 

One of the main objectives of the investigation was to ascertain 
whether motions through the stall and into the incipient spin could be 
reliably indicated by means of free dynamic models. The repeat test, 
the results of which are shown in figure 5, was made with the test con- 
ditions nearly the same, insofar as the available equipment allowed, 
as those of the original or initial test in order to determine whether 
the motion would repeat itself. Although there were some differences 
for the repeat test in that the model had a more positive pitching 
velocity and a somewhat different sideslip angle than for the initial 
test, the results in general are in good agreement. Quantitatively, 
the results are generally close to or within the estimated limits of 
accuracy and indicate that the model probably repeats its motions for 
similar launching conditions. Because of the unsteady nature of the 
air flow which exists at and beyond the stall, it has been thought that 
airplanes may not repeat motions in the stalled-flight regime. The 
results from the repeat test and experience in spin and spin-recovery 
work, however, indicate that motions at and beyond the stall may gen- 
erally repeat themselves except for critical cases such as have been 
encountered in studies of developed spins when the spin and recovery 
change unaccountably. 

The motions discussed are for power-off stalls and the use of power 
may appreciably alter the motion or the symmetry from that presented 
herein. The results should not be interpreted as representing the 
motions of an airplane corresponding precisely to the model, because of 
scale effects, but can be considered indicative of the trends in the 
motion of a typical personal-owner or liaison airplane. Changes in 
configuration, wing loading, and mass distribution also may alter the 
results appreciably. 

In order to evaluate the symmetry of the model and to ascertain 
that the motion to the right previously discussed was not influenced 
primarily by model asymmetry, tests were conducted to determine the 
motion with the rudder set full left ( 3 O 0 ) • The results are presented 
in figure 6 . As before, the variation of angle of attack, angle of 
sideslip, and angular velocities with time are presented as the model 
stalled and rolled to the left prior to entering a spin to the left. 
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A comparison of figure 6 with figure 5 indicates that the model, which 
appeared to be dimensionally symmetrical, showed generally similar motions 
when the rudder was set either to the right or to the left. 


CONCLUDING REMARKS 


An investigation has been conducted to study the motion of a free 
dynamic model of a personal-owner or liaison airplane through the stall 
and into the incipient spin. This dynamic-model technique does not lend 
itself to specific or development studies of a given design because of 
probable scale effects existing at the stall; however, the method is 
applicable for indicating the trends of the motion through the stall and 
into the incipient spin and possibly the effect of certain parameters on 
this motion. The results showed that the model became unstalled in 
inverted flight after its initial rotation of substantially l 80 ° in roll. 
Since the rates of rotation were low, it may be inferred that termination 
of the incipient spin would be most readily obtained, by proper control 
movements at this point. The results of the model flights show that 
the motion is generally repeatable and that, power effects being excluded, 
a geometrically symmetrical model should have fairly symmetrical motions 
to the right and to the left. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., January 28, 1953- 
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APPENDIX 

DERIVATION OF EQUATIONS FOR EVALUATION OF 
PARAMETERS IN INCIPIENT-SPIN MOTION 


The determination of angle of attack, angle of sideslip, and the 
angular velocities of the model in flight are based upon calculating 
its component velocities along the body axes. The determination of 
these velocities is accomplished by differentiating the respective 
earth-axis space coordinates of the model against time to obtain the 
earth-axes velocities of the model and then resolving by trigonometric 
means these velocities to a body-axes system. The following process 
was used. 

If X, Y, and Z are the earth axes, and Xg, Yg, and Zg are 
the body axes, both systems having their origin at the center of gravity 
of the model, the relationships between earth-axes velocities and body- 
axes velocities are as follows: 


ug = u cos {_ XXg + v cos [_ YXg + w cos /_ ZXg (Al) 


Vg = u cos J_ XYg + v cos [_ YYg + w cos j_ ZYg (A2) 


wg = u cos j_ XZg + v cos J_ YZg + w cos [_ ZZg (A 3 ) 


where [_ XXg, J_ YXg, and so forth, are angles between the earth and 

body axes designated. These nine angles, which relate the velocities 
along the earth and body axes through direction cosines, therefore must be 
found. The coordinates of the model right wing and tail reference points 
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(wing tip in Xg,Yg plane and tail cone in Xg,Zg plane, respectively) 
from the center of gravity are determined: 


x R = x_ - x_ 

^ ~R T cg 


XR X r ~ Xp 

R eg 


Z'D — Zy. “ Z-m 

K r R r cg 


= - ^cg 


- Xr T - yr cg 


2T - Zr T " 2r cg 


From the foregoing coordinates, the angles 0, 0', 0 ' 

and 0, as well as i|r", can be obtained (see fig. l) : 


z T 

0 = sin" 1 — 

Z<ji 

y T 

i|r = tan” 1 — 

Xrp 

0 ' = sin -1 -2- 
b/2 

i|r 1 = tan -1 — — 

XR 


V = ♦’ - * 


0' = cos _1 (cos 0 cos 1 ^") 




(A4) 


0 


cos _1 (cos 0' cos \jr") 
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From figures 1 and 8 and from the determination of the angles cal- 
culated in equations (A4 ) , the trigonometric functions necessary for the 
solution of equations (Al) and (A2) can be obtained as follows: 

cos l_ XXjj = cos \|r cos 0 (A5) 

cos j_ YX-g = cos(90 - \|r)cos 9 = sin \|r cos 0 (A6) 

cos / ZXb = cos (90 +0) = -sin 0 (A7) 

cos j_ YYg = cos i|r' cos 0' (a8) 

cos j_ XYg = cos(90 + \lr')cos <f' = -sin \(r ' cos 0' (A9) 

cos j_ ZYb = cos (90 - 0') = sin 0' (AlO) 

The determination of the trigonometric functions needed for the 
solution of equation (A 3 ) is somewhat more difficult than those for 
equations (Al) and (A2) . These functions may be found by consideration 
of the spherical triangles depicted in figure 8. 

cos ]_ ZZg = cos 0 cos 0 = cos 0' cos 0' (All) 

Also, from figure 8 

cos J_ XZg = sin [_ ZZg cos [_ XZZq 
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where 


sm 


L ZZ B = \fl - 


cos z 0 


COS 


20 


and 


cos L XZZ B = cos \Jr cos l_ ZgZF + sin \|r sin Z^ZF 


or 


cos J_ XZZ-g = 


cos \|r sin 9 cos 0 + sin \|r sin 0 

^1 - cos^0 cos^0 


Therefore, 


cos £ XZrg = cos ^ sin 6 cos 0 + sin \|r sin 0 


(A12) 


Now 


cos [_ YZ3 = sin l_ ZZ3 cos J__ YZZg 


where 


cos Z yzz b = sin \|r cos j_ Z^ZF - cos \|r sin /_ ZrgZF 


or 


cos / YZZ b = 


sin \Jr sin 0 cos 0 - cos sin 0 
\j 1 - cos ^9 cos^0 
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Therefore , 


cos l_ YZ-g = sin \|r sin 0 cos 0 - cos \|r sin 0 (AI 3 ) 

(Note that angles XZZg, YZZg, and ZgZF are angles between the arcs 

of the spherical triangles shown in figure 8.) Thus, with the equa- 
tions (A5) to (AI 3 ), the nine necessary direction cosines are determined 
and equations (Al) to (A 3 ) may be solved for the velocities along the 
body axes . The velocities along the earth axes necessary for the solu- 
tion of these equations are obtained, as was indicated previously, by 
taking the first derivative of the space coordinates with respect to 
time. 


From a solution of equations (Al) and (A 3 ) wherein a determination 
of the velocities of the center of gravity along the X and Z body 
axes has been made, the angle of attack, measured at the center of gravity 
and defined as the angle between the X body axis and the projection of 
the relative wind on the Xg,Zg plane, can be calculated as follows: 


a 


w-o 

= tan~l — 
U B 


(A14) 


If the resultant velocity of the model is calculated 


V = 



+ v- 


B 


+ w- 


B 


(A15) 


the angle of sideslip may be calculated as follows: 


P 


sin“l 


2b 

V 


(A16) 


For a determination of the angular velocities p, q, and r, refer- 
ence 4 was used. The angular velocities are determined on the basis of 
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Euler's angles, which are herein defined as \|r, 9 , and 0. The equa- 

tions for the angular velocities therefore are as follows: 


p = 0 - sin 9 


(A17) 


q =9 cos 0 + ^ sin 0 cos 9 


(A18) 


r = ijr cos 0 cos 9-9 sin 0 


(A19) 
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TABLE I— DIMENSIONAL AND MASS CHARACTERISTICS OF THE 
CORRESPONDING FULL-SCALE AIRPLANE 


Dimens ional: 

Over-all length, ft 22.37 

Wing: 

Airfoil section MCA 43012 

Incidence, deg 3 

Dihedral, deg 6 

Twist, deg 0 

Span, ft 33*^3 

Mean aerodynamic chord, c, ft 4.89 

Leading edge of mean aerodynamic chord rearward of leading edge of wing, ft 0.05 

Taper ratio 1.00 

Area, sq ft 163.28 

Aspect ratio 6.93 

Ailerons: 

Span, ft 7.15 

Area, rearward of hinge line, sq ft 15.70 

Horizontal-tail surface: 

Span, ft 10.25 

Total area, sq ft 26.39 

Elevator area rearward of hinge line, sq ft 11.02 

Aspect ratio 3*98 

Incidence, deg 0 

Dihedral, deg _ 0 

Distance from quarter chord of c to elevator hinge line, ft 13*73 

Section Modified NACA 0009 

Vertical-tail surface: 

Span, ft 5 *32 

Total area, sq ft 12.96 

Rudder area rearward of hinge line, sq ft 6.48 

Aspect ratio 1.26 

Offset, deg 0 

Distance from quarter chord c to rudder hinge line, ft 14.18 

Section Modified MCA 0009 

Mass: 

Weight, lb 2449 

Wing loading, lb/sq ft 15*0 

Relative density: 

Sea level 5.82 

5,000 ft 6.76 

Center of gravity: 

x/c 0.26 

z/l 0-12 

Moments of inertia (about center of gravity) , slug-ft 2 

% 

Iy B 1998 

^ 

Inertia parameters: 

IX B ' lY B 45 

mb 2 

lY B - IZ B -251 

mb 2 

IZ B ' IX B 206 

nib 2 ^NACA: 
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Figure 1.- Diagram of earth and body axes through the center of gravity 
of the model showing angular relations . Positive directions of axes 
and positive values of angles are shown except for 0 and 9 ' which 
are shown negative . 
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Figure 2.- Three-view drawing of model tested. 
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Figure 3 .- Photograph of model tested. 




Figure k.- Photograph of launching apparatus. L-76236.I 
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(a) Angle of attack, angle of sideslip, and altitude loss. 

Figure 5.- Time histories of some components of model motion for the 
original and repeat right roll-offs. 
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(b) Rolling, pitching, and yawing velocities. 


Figure 5.- Concluded. 


Angle of attack, a, deg Angle of sideslip, $ , deg Altitude loss, h, ft 


2k 
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(a) Angle of attack, angle of sideslip, and altitude loss. 

Figure 6.- Time histories of some components of model motion for the 

symmetry-check left roll-off. 


Pitching velocity , q , radians /sec Yawing velocity , r , rodions/sec 
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(b) Rolling, pitching, and yawing velocities. 
Figure 6.- Concluded. 
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Figure 7 .- Concluded. 
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Figure 8.- Diagram of spherical triangles used in the analysis of data. 
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